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Abstract 



Lattice measurements of the Pisa group (A.Di Giacomo et al., hep-lat/9603018 ) 
psj \ are analyzed numerically and parameters of correlation functions are extracted from 

£H ■ the data - both below and above deconfinement temperature T c . Gluon condensate is 

found for six temperatures in the interval 0.956 T c - 1.131 T c and field distributions in 
deconfined phase are obtained. 

<n : 

: 1 Introduction 

> . 
o . 

. Study of field distributions around static quarks has a long history. The information obtained 
both analytically and on the lattice has an important meaning in several respects. Firstly, it 
\ demonstrates clearly the appearance of the QCD string between the charges in the confining 
phase and the detailed contents of the fields in this string; e.g., the string consists mainly 
of longitudinal color-electric field. Secondly, since the string enters in many dynamical 
quantities, such as static interquark potential, spin-dependent forces for heavy quarkonia 
etc., one can easily compute these quantities from the string field distributions. 

Finally, and this is the purpose of the present paper, the phenomenon of deconfinement 
^ ' is not fully understood from the point of view of the string and field distributions. 

What actually happens when temperature T exceeds the critical value - the string disap- 
^ ■ pears or distribution of fields drastically change, so that to compensate the string tension? 

Happily there exist numerical lattice measurements of field correlators near the critical 
temperature T c , made by the Pisa group [1], where both electric and magnetic correlators 
are found with good accuracy. These data clearly demonstrate the strong suppression of 
color-electric component above T c and persistence of color-magnetic components. 

The purpose of the present paper is twofold. First, we reanalyze the data of [p]] in 
terms of correlators D ■ and D l ' which are better understood from the point of view of 
perturbative and nonperturbative contributions [2]. There was e.g. shown that the string 
tension a ("conventional" electric or "spatial" magnetic) is expressed directly through the 
D E ' B but not the D E ' B . Given the simple assumption of behavior D E ' B (x) and D E,B (x) in all 
region of distances x between correlation points, we obtain the behavior of gluon condensate 
near the T c (which is defined by D E ' B (0) and D E,B (0)). Second, using the obtained D E,B {x) 
and D E ' B (x), we calculate the color field distributions around static quarks in the deconfined 
phase. We perform these calculations using the connected probe [3] in the framework of Field 
Correlator Method (FCM) [4,5]. We study in detail two possible from the deconfinement 
point of view regimes, corresponding to two forms of D(x), extracted from lattice data [1]. 
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In the first one the string disappears and in the second (less appealing physically, but more 
supported by the lattice data) the string becomes a coaxial cable with the empty core. 
The paper is organized as follows. 

In Section 2 a short description is given of our fitting procedure of magnetic correlators 
at any T and the electric correlators at T < T c , while in Section 3 more prolonged one is 
given for electric correlators at T > T c . The Section 4 is devoted to the behaviour of the 
gluonic condensate below and above T c . In Section 5 the detailed equations for the field 
distributions around deconfined quark and antiquark are given and illustrated graphically. 
In the concluding section results of the paper are summarized and discussed. 

2 Fitting data with nonzero string tension 

In this section we consider magnetic correlators in all temperature region and electric cor- 
relators at T < T c since both of them produce (e.g., spatial) nonzero string tension. We fit 
the data [1] using the method of least squares [6]. To begin with, we express the functions 
D\\(x), D±(x) through the D(x), Di(x) and represent the lasts as sums of nonperturbative 
(NP) and perturbative (P) (diverging at zero) contributions. This is the procedure used by 
A.Di Giacomo et al.[7] when analyzing the correlation functions at zero temperature. The 
difference is that we should distinguish electric and magnetic correlators due to the fact that 
the finite temperature theory has only the 0(3) x 0(1) symmetry. 
We fit distributions of 

D^(x) = D*>»(x) + Df*(x) + X -^j^ (1) 

and 

D^ B (x) =D E > B (x) + D*' B {x) (2) 

electric and magnetic correlation functions in the range from 0.4 to 1 fm. All points are 
measured at x 4 = 0. In this section we parameterize the functions as follows: 

D(x) = Aexp(-x/X A ) + -^ z exp(-x/X a ), (3) 

b 

D^x) = Bexp(-x/\ B ) + — exp(-x/\ b ). (4) 

ar 

Fitting the magnetic functions (which are close to exponentials), we set in (3), (4) A a = 
A& = As, thus using k = 6 fitting parameters. Given N = 12 number of data, we have as a 
result n = N — k = 6 number of degrees of freedom. The results are shown in Table 1 . One- 
standard- deviation errors Actj are determined from the equation x 2 ( a ii ■■■■> «i + Aa:j, a^) = 
X 2 («i, ccj, a k )+Axl, where are fitted parameters and Axf = 1, Axl — 4.75, Axl = 7 
[6]. 

Preliminary fitting of electric data at T = 0.956T C with parameters of (3), (4) gave un- 
acceptable big x 2 due to the end points of and Xa, A a , and \ b were found close to 
each other. Therefore we had removed mentioned points and set A^ = A a = A& to get the 
reasonable fit (Table 2). 

Preliminary fitting of electric data at T = 0.978T C have shown that x 2 is a l so too big 
and besides that A^ is close to Aj, and A# to A a . Therefore we set Aj, = A^ and A a = A#- 
To improve % 2 , we enlarge two times the error of the last point of (Fig. 1). One can 
make sure from Figure 1 that this point is largely off the exponential curve, which may be 
connected to the lattice size effects. The results of fitting are shown in Table 2. 
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3 Fitting electric data in the deconfinement region 



above T c has a drop that is presumably related with the deconfinement transition, when 
the string tension of area law asymptotics of Wilson loop for static quarks disappears. In 
gaussian approximation of FCM [5] we get 

a = i J dxidxiD^ { Xl , x 4 ) = 0. (5) 

In this region we should alter the form of D (3) to justify (5). It is naturally to set 

D NP (x) = 0. (6) 

(We omit subscript "E" here and in what follows in this section.) As we shall see below, this 
form assures reasonably well fit of data. Let us propose another form of D to ensure good 
fitting of data. We suppose now more constrained 0(4) symmetry: -D(x 2 , xf) = D(y?+xf) = 
D(x 2 ) to get 

a = | dx 2 D^{x 2 ) = | dx 2 Df{x 2 ) = 0. (7) 
To better reproduce the mentioned drop of data, we set D^ p (x) = 0, i.e., 

£NP {x) = _ d np {x) _ x dDf p (x) ^ (g) 



2 dx 



leaving meanwhile D\ in form (4) intact. 
So far we have adopted two forms of D(x) 



D(x) = -^exp(-x/X a ); (9) 



x 4 



D(x) = - - 1) exp(-x/A B ) - 4 exp(-x/A ). (10) 

2Ab 



We fit data on Du and D± at T = 1.011T C in two ways, using (a): functions (9), (4) and (b): 
(10), (4). Having in mind the small number of data we set A a = A& = A# to have k — 4. 
Given = 9, we obtain n = 5 (Table 3). We see from the Table 3 the somewhat reasonable 
X 2 /n = 1.7 in case (a) and excellent x 2 / n — 1-05 in case (b). 

At higher temperatures there are only two measurements of D\\, with values significantly 
less even than the errors of corresponding points of D±. This circumstance allows us to 
subdivide the fitting procedure in two stages. At the first stage we fit the difference D± — D\\ 
(cf. (1),(2)) by —x 2 dD 1 /dx 2 (cf. (4)), and extract parameteres (B, \s,b, A&). At the second 
stage we fit D±, reproducing it (in two ways, due to two cases of D) by the sum of D and Di 
with parameters extracted from the first fit; D is taken in forms (9) and (10) with A a = A& 
in both cases, and only parameter a is allowed to vary (Table 4, Fig. 2). One-standard- 
deviation errors are determined as described above, with k = 4 at the first stage and k — 1 
at the second stage of fitting. In the Table 4 the second-stage-fitting results are separated 
by horizontal line; x\i refer to the first and second stages correspondingly. 

As Dj_ — D\\ Pi D±, at the first stage we actually fit the data by 

Bx 26 x 

D±(x) = ^- ex P (-x/\ B ) + _(! + _) ex P (-x/\ b ) (11) 
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and at the second stage in case (a) by 

D ± (x) = B exp(-x/X B ) + exp(-x/\ b ) (12) 

Ob 

and in case (b) by 

Bx a ~\- b 

D±(x) = — — exp(-x/\ B ) H j- exp(-x/\ b ), (13) 

1\b x 

with all parameters of (12), (13) except a fixed by the first stage. One could see that (13) 
well reproduces (11) at a = b and any B and Ag, for i <C 4A&, i.e., in all measured region 
0.4 fm< x < 1 fm (see Table 4). 



4 Temperature dependence of the gluon condensate 

The gluon condensate is defined as 

G2 = t(W, (14) 



7T 



where a s is the strong coupling constant, are the gauge field strengths taken at the point 
x = and the averaging is performed over all vacuum configurations. At zero temperature 
the FCM reads 



92 (C(^(n^?(^,^)^w(^)^i(^^)> = 

2\ , n tu2\ 



(V- - Spv^)(D(h 2 ) + zm/^H 



dDi(h 2 ) 

(h^h p 5 ua - h^h a 5 up - h p h u 5^ + h u h a 5 pp ) — — 2 — , (15) 

where h = x — x'; a,. ..,5 are color indices. At x = x' = one uses in (15) tiT a T b = \8 ab 
and gets 

<j2 (F a F a ) = 12(D(0) + D 1 (0)). (16) 



In what follows we shall use the zero temperature lattice results [7]: 

£> NP (0) = 3.3 x 10 8 A£ = 129 fm" 4 , Df p (0) = 0.7 x 10 8 A 4 = 27 fm" 4 , (17) 

where A^ = 0.025 fm -1 is the fundamental constant of QCD in lattice renormalization 
scheme; its value is extracted from the string tension. 
At finite temperature one derives from (15) 

//2 (^(0)^(0)) = 3(D E (0)+Df(0)), (18) 



2N, 
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;? _ Vc (5?(0)5?(0)> = 3( J D B (0) + J D B (0)), (19) 

where Ei = F 0i and Bi = \eu m Fi m . Note that at finite temperature -D E and _D B acquire 
subscripts for symmetry breaking 0(4) — > 0(3) x 0(1) reason. We have to distinguish 
electric and magnetic contributions to condensate: 

(F; v F; v ) = (F 2 ) cl + (F 2 ) magn , (20) 



where 

(F 2 ) cl = (F«F«) + (F»F«) = 2{E?E?), (21) 

(i 7l2 )ma g n = -< iji.< ij I, {F, h Fjh) + -^ ijk e ijk {F k:j F kj ) = 2{B°-B°-) (22) 

Substituting (16), (18)-(22) into (14), one obtains that normalized gluonic condensate is 

G 2 (T) _ D E - NP (0) + D[ ' NP (0) + D B - NP (0) + D* ' NP (0) 
G2 " 2(D(0)+D 1 (0)) ' ( } 

where G 2 (T) is the gluon condensate at temperature T; G\ = (^(O). 

According to our fitting, magnetic condensate is D B ' NP (0) ' NP (0) = A + B, with A, B 
taken from Table 1. Electric condensate at temperature T < T c is D E ' NP (0) + -Df' NP (0) = 
A + B, with A, B taken from Table 2. At T > T c in case (a) £> E ' NP (0) = 0, L>f' NP (0) = B, 
with 5 taken from Tables 3,4, and in case (b) D E > NP (0) + D? ,NP (0) = 0. 

The behavior of the condensate and that of its electric and magnetic constituents with 
temperature are shown in Tables 5,6. Data on the whole condensate from these Tables are 
plotted in Fig. 3. We see that at T < T c the value of the condensate is close to its zero 
temperature value. At T > T c in the case (a) there is a fast growth of condensate. We will 
discuss its physical meaning in the concluding section. In the case (b) the condensate value 
is about half of its zero temperature value. 



5 Field distributions around deconfined quarks 

In this section we consider NP part of gluodynamical field generated by static QQ sources, 
using the connected probe [3,8]. There was shown that the only nonzero components in 
this system are longitudinal (along quark axis) and transverse electric fields Ei(xi,x 2 ) and 
E 2 (xi,X2), where x\ is coordinate along quark axis and x 2 - distance to the axis. 

In case (a), when D NP = 0, NP part of D\ is D\ = B exp(— x/X), where x = \J x\ + x\ + x\ 
(due to axial symmetry we may set x% = 0), A means Ag. Here and in what follows we omit 
subscript " E" . From the equations of [8] 



(E 1 (x 1 ,x 2 )} i ^ Q = J* dx[ J°^dx' 4 (D^h 2 ) + (hi + hi) 



dD 1 (h 2 
dh 2 



Bix^xl + xlK^xj + xl/X) - ( Xl - R)y/(x 1 - R) 2 + x 2 ^ { Xl - Rf + x 2 /X)), (24) 

rR ,oo dD^h 2 ) 



tit roo 

(E 2 (x 1 ,x 2 ))qq = J o dx[J dx' A hih 2 



dh 2 

Bx^xl + xlK^xl + xl/X) - y/fa - Rf + x 2 2 K 1 { y /{x 1 - R) 2 + x 2 /X)), (25) 
where K\ is McDonald function. The total field is 

(E) 2 = (E,) 2 + (E 2 ) 2 . (26) 

In Fig. 4 we plot (E(xi, x 2 )) 2 distribution with parameters corresponding to the case T = 
1.070T C and R = 2 fm. We observe two "volcanoes" with quarks hidden in their bottoms. 
These two spherically symmetrical in coordinate space distributions are defined as 

E{r)=Br 2 K 1 (r/\), (27) 
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where r is a distance from quark or antiquark. The field at quark and antiquark positions 
is zero and linearly rises at small r. Maximal value of field is 



£™ = £(1.33A) = 0.63BA 2 . (28) 

In Fig. 5 the vector field distribution (E(xi,X2)) is shown in the vicinity of the quark. 
In the case (b) the transverse part of the field, E 2 , remains the same (25). Let us calculate 
using (8) the longitudinal part of the field, E\\ 

m rR r°° ( f)D-,(h 2 



rR /-oo / 

/ dx[ / dx' A D(h 2 ) + D^h 2 ) + {h\ + h\) 

JO J~oo \ 



dh 2 



dhi / dhi(-h 2 )—— = dhi dh^—Dx 

xx-R J-oo Oh z J Xl -R JO h\ 



Bx 2 r Xi 



dh.K^hl + xl/X). (29) 



A Jxx-R 

In Fig. 6 we plot (29) for T = 1.070T C and R = 2 fm to observe the "double quasistring" . 
The quasistring profile, Ei(x 2 ) = (Ei(R/2, x 2 ))qL at R — > oo, is 

Ex{x 2 ) =7rSa^exp(-|x 2 |/A). (30) 
The field in the centre of quasistring is absent. The maximal value of the field is 

^max = El {2\) = 1.7BX 2 . (31) 

In the coordinate space the quasistring resembles coaxial cable with empty core and tube 
shell. In Fig. 7 we plot (E(xi, a^))^ distribution around Q and Q. 



6 Conclusions 

Results of our paper based on the analysis of the lattice data on correlation functions at 
finite T [1] give a full support of the dynamical picture of deconfinement, which was first 
suggested in [9]. 

Namely, confining and deconfining phases according to [9] differ first of all in the vacuum 
fields, i.e., in the value of the condensate and in the field correlators. It was argued in [9] that 
color magnetic correlators and their contribution to the condensate are kept intact across the 
temperature phase transition, while the confining electric part abruptly disappears above T c . 
Both features are present in [1] and in the results of the present paper. Indeed, one can see 
from Tables 5,6, that the magnetic part of condensate is roughly constant around T c . The 
role of magnetic field above T c was mentioned repeatedly in the literature, see recent lattice 
reviews [10], it reveals itself in particular in creating nonzero spatial string tension and 
so-called screening hadronic masses - (see |TT| and refs. therein). 



The situation with electric fields is more subtle, as can be seen from our results. In 
two possible situations have been considered when nonperturbative part of D E vanishes or 
stays nonzero above T c . From Figs. 4,6 one can clearly see the field distributions in cases 
of two possible solutions, (a) and (b), one with vanishing _D E,NP , another with nonvanishing 
but oscillating _D E,NP , both yielding zero string tension in the deconfinement region. 
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In the ) the electric contribution to condensate is determined by D 1 ' (0). While 

_D E ' NP (x) vanishes identically, the Z}f' NP correlator is different from zero above T c and its 
contribution to the condensate grows sharply with temperature. Hence the role of Z}f' NP 
appears in creating the sharp rise of the " nonideality" (e — 3p) just above T c (cf. Fig. 3). One 
special remark is due to the regime (b), where D E,NP (x) is nonzero above T c but changing 
sign. This regime creates rather peculiar picture of fields — the " quasistring" with the empty 
core and surrounding it tube shell at two correlation length distance from the quark axis 
(Fig. 6). 

As the regime (a) seems to be more natural from physical point of view, one should study 
in more detail the consequences of the strongly increasing with T Z)f ,NP (0) in the deconfining 
region. 

The impossibility of resolving our present ambiguity (regimes (a) and (b)) calls for fur- 
ther numerical and analytical studies. It is necessary for understanding of the dynamics 
of the phase transition, where Polyakov loops and hence color-electric fields may play very 
important role. 

The authors are grateful to A.Di Giacomo for useful remarks and suggestions; the partial 
support of grants 00-02-17836 and 00-15-96786 is gratefully acknowledged. 
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List of tables 



Table 1: Parameters of D B . 



T 


0.956T C 


0.978T C 


1.011T C 


A, fm~ 4 


188.9±2.4 


154±2 


183.1±3.0 


X A , fm 


0.1917±0.0007 


0.2047±0.0008 


0.1852±0.0008 


B, fm~ 4 


7.7±0.6 


7.7±0.8 


1.87±0.13 


Xb, fm 


0.380±0.008 


0.344±0.007 


1.11±0.04 


a 


1.11±0.06 


1.46±0.07 


0.64±0.03 


D 


n q A _i_n no 
0.34±0.03 


n A ^7_i_n r\A 
U.4w±0.04 


n Q£_i_n no 
0.3o±0.Uz 


X 2 / n 


0.62 


1.28 


1.43 


T 


1.034T C 


1.070T C 


i.i3ir c 


A, fm~ 4 


128.8±2.3 


111.9±2.0 


150.9±2.5 


Xa, fm 


0.210±0.001 


0.2191±0.0012 


0.2009±0.0010 


B, fm" 4 


2.92±0.23 


3.54±0.29 


3.07±0.23 


A B , fm 


0.69±0.02 


0.631±0.017 


0.774±0.024 


a 


0.92±0.04 


1.039±0.043 


0.885±0.031 


b 


0.47±0.03 


0.525±0.032 


0.506±0.023 


X 2 / n 


0.53 


1.18 


0.58 



Table 2: Parameters of D E below T c . 



T 


0.956T C 


0.978T C 


A, fm" 4 


228±4 


189.5±4.2 


Xa, fm 


0.1823±0.0008 


0.1812±0.0011 


B, fm" 4 


10.8±0.5 


14.3±0.6 


Xb, fm 


0.435±0.011 


0.411±0.006 


a 


3.1±0.4 


0.99±0.08 


b 


0.9±0.2 


0.71±0.18 


X 2 / n 


0.39 


1.7 
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Table 3: Parameters of D E at T = 1.011T C . 





(a): D N,J = 


(b): Df^ = 


B, fm- 4 


2.46±0.24 


1.96±0.34 


As, fm 


0.682±0.016 


1.072±0.047 


a 


1.74±0.05 


1.36±0.04 


b 


0.86±0.04 


0.68±0.03 


X 2 / n 


1.7 


1.05 



Table 4: Parameters of D E above T c . 



T 


1.034T C 


1.070T C 


1.131T C 


B, fm~ 4 


80±19 


235±53 


519±99 


X B , fm 


0.121±0.05 


0.105±0.004 


0.37±0.03 


b 


0.63±0.03 


0.48±0.04 


0.41±0.03 


Aft, fm 


0.86±0.05 


1.0±0.1 


0.86±0.08 




0.69 


0.9 


0.028 


(a) a 


1.21±0.04 


0.94±0.04 


0.86±0.03 


(a) x| 


1.26 


0.9 


2.9 


(b) a 


0.96±0.04 


0.66±0.04 


0.61±0.03 


(b) xl 


1.66 


0.23 


0.56 



Table 5: Gluonic condensate below T c . 



T 


0.956T C 


0.978T C 


L> B > NP (0) + Df > WF (0), fm" 4 
L»e,np (0) + ^e,np (0)) fm _4 

G 2 (T)/Gl 


196.6±2.5 
238.8±4.0 

1.393±0.015 


161.7±2.2 
203.8±4.2 

1.171±0.015 



Table 6: Gluonic condensate above T c . 



T 


1.011T C 


1.034T C 


1.070T C 


1.131T C 


D B ^ p (0) + D*' Ntl (0), fm- 4 

G 2 (T)/G° 2 (a) 

(b) 


185±3 

0.60±0.01 
0.59±0.01 


131.7±2.3 

0.69±0.06 
0.422±0.007 


115.4±2.0 

1.12±0.17 
0.370±0.006 


154.0±2.5 

2.16±0.32 
0.493±0.008 
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X 



0.5 0.6 0.7 0.8 0.9 1 1.1 X 



Figure 1: Fitted functions h\D^(x) and \nDf(x) (shown by solid lines) at T = 0.978T C . x 
is measured in fm. and D\ are measured in fm -4 . Data are shown by points. Errors 
of data are comparable with size of the points. The separately shown error is two times 
enlarged by hand for the improvement of the fit. 
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In D E , 




x 



Figure 2: Fitted function ha.D^_(x), shown by solid line in case (a) and dashed line in case 
(b) at T = 1.070T C . x is measured in fm. is measured in fm~ 4 . Data are shown by points 
joined by dashed lines. Errors of data are comparable with size of the points. 



G 2 /G^ 
2.5 




0.975 1.025 1.075 1.125 



Figure 3: Gluon condensate dependence on T/T c near the deconfmement transition. Con- 
densate is measured in its zero temperature value units. In the deconfmement region the 
magnetic part of condensate does not considerably change. Electric fields give rapidly rising 
contribution to condensate in case (a) and do no contribute to condensate in case (b). 
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Figure 4: The total field distribution (E(xi, x 2 )) 2 , measured in fm 4 , in case (a) at T = 
1.070T C . x\ and £2 are measured in fm. QQ separation R = 2 fm. 




Figure 5: Vector distribution of (E(xi,X2)) in case (a). —0.5 fm< x\ < 0.5 fm, —0.5 
fm< x 2 <0.5 fm. Quark position is marked with disk. Points of maximal value of field are 
marked with dashed circle of radius 1.33A. 
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Figure 6: Distribution of {E 1 (xi,x 2 )), measured in fm 2 , in case (b) at T = 1.070T C . 



^ +~ - »- *~ — *~ - - fc. - ^ ^ ^ 




Figure 7: Vector distribution of (Ef^i,^)) in case (b). Positions of Q and Q are marked 
with disks. —0.1 fm< X\ <2.1 fm, —0.5 fm< x 2 <0.5 fm. 
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